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MOLECULAR  MODELING:  AN  APPROACH  FOR 
THE  STUDY  OF  PIEZOELECTRIC  POLYMERS 


INTRODUCTION 

The  central  theme  of  this  report  is  to  demonstrate  how  molecular  modeling 
can  be  used  to  study,  predict,  and  tailor  the  piezoelectric  properties  of  semicrystalline 
polyvinylidene  fluoride  (PVDF)  and  its  copolymer  with  trifluoroethylene  P(VDF- 
TrFE).  This  is  an  application  where  molecular  modeling  has  been  successfully  used 
to  understand  and  improve  the  properties  of  piezoelectric  polymers  for  active  and 
passive  sonar.  PVDF  and  P(VDF-TrFE)  are  remarkably  versatile  polymers.  PVDF 
consists  of  thousands  of  repeating  monomeric  units  of  vinylidene  fluoride  (VDF). 

Since  this  polymer  is  prepared  from  a  single  monomer,  it  is  often  referred  to  as  a 
homopolymer.  A  schematic  of  PVDF  is  shown  in  figure  1  where  "A"  represents  the 
VDF  monomer.  P(VDF-TrFE)  is  composed  of  vinylidene  fluoride  and 
trifluoroethylene  (TrFE)  monomers.  These  monomers  can  be  arranged  in  varying 
ways  leading  to  several  possible  structures.  P(VDF-TrFE)  is  said  to  be  a  random 
copolymer  since  the  distribution  of  its  monomers  is  random.  A  possible  arrangement 
is  shown  in  figure  1  for  "A"  and  "B"  monomers.  Since  copolymer  compositions  may 
differ  in  their  fluorine  content,  the  properties  of  these  compositions  are  determined 
by  the  relative  ratios  of  the  VDF  and  TrFE  monomers. 

Both  PVDF  and  P(VDF-TrFE)  can  be  processed  into  semicrystalline,  nonpolar 
materials.  When  poled  in  a  strong  electric  field,  P(VDF-TrFE)  is  converted  into  a 
polar,  piezoelectric  material.^  PVDF,  however,  must  be  stretched  before  it  converts 
into  its  piezoelectric  form.  When  mechanically  stressed,  piezoelectric  PVDF  and 
P(VDF-TrFE)  generate  an  electric  polarization;  when  electrically  stressed,  they 
undergo  mechanical  deformation.^”^  The  polarization  of  PVDF  and  P(VFE-TrFE)  is 
due  to  a  highly  oriented  alignment  of  the  numerous  carbon-fluorine  dipoles  along  the 
direction  of  the  electric  field.  The  position  of  fluorine  in  the  periodic  table  guarantees 
its  strong  electronegativity,  which  means  that  it  can  attract  electrons  from  the 
neighboring  carbon  atom  to  which  it  is  attached.  This  leads  to  the  development  of  a 
carbon-fluorine  dipole.  The  presence  of  fluorine  also  confers  PVDF  and  P(VDF-TrFE) 
polymers  with  their  remarkable  properties  including,  for  example,  their  hydrophobic- 
ity,  water  insolubility,  stability  to  thermal  and  oxidative  degradation,  and  resistance 
to  fire,  oil,  and  organic  solvents. 

Piezoelectric  fluoropolymers  offer  several  advantages  when  compared  with 
traditional  piezoelectric  ceramics  used  in  underwater  sonar.  Polymers  are  light, 
mechanically  flexible,  robust,  conformal,  and  better  acoustically-matched  with  water 
than  conventional  ceramics.^  Because  of  these  features,  PVDF  has  been  considered  as 
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a  passive  hydrophone  material  for  large  aperture  arrays  on  future  submarines. 
P(VDF-TrFE)  is  a  potential  candidate  for  underwater  acoustic  projectors  of  broad 
bandwidth  requiring  stacks  of  thick  electroded  elements.^  Since  there  is  no  require¬ 
ment  for  mechanical  stretching  of  the  copolymer  chains,  solvent-cast  or  melted 
copolymer  can  be  poured  into  molds,  annealed  to  induce  crystalli2ation,  and  then 
poled. 


At  the  outset  of  this  project,  the  question  of  polymer  composition  emerges  as 
an  issue  of  central  importance.  For  example,  how  do  the  mechanical,  electrical,  and 
thermal  properties  vary  as  a  function  of  fluorine  content?  In  the  case  of  P(VDF- 
TrFE),  what  are  the  optimum  compositions  for  piezoelectric  performance?  Can  a 
predictive  model  be  formulated  to  allow  the  transducer  designer  to  select  the  most 
appropriate  composition  for  a  specific  sonar  application?  The  objective  of  this 
research  was  to  use  molecular  modeling  to  gain  deeper  insight  into  the  piezoelectric 
properties  of  PVDF  and  P(VDF-TrFE)  and  to  develop  a  practical  predictive  frame¬ 
work  to  assist  the  transducer  designer  in  selecting  and  tailoring  the  properties  of 
various  polymer  compositions  for  underwater  sonar  transduction  applications. 


THEORY 

THE  MOLECULAR  AND  CRYSTAL  STRUCTURES  OF  PVDF  AND 
P(VDF-TrFE) 

In  this  section,  we  describe  briefly  the  structural  characteristics  of  PVDF  and 
P(VDF-TrFE).  PVDF  is  synthesized  by  the  polymerization  of  1,1-difluoroethylene,  or 
vinylidene  fluoride.^  The  repeating  monomer  units  have  the  formula  -CH^-CFj-, 
where  C,  H,  and  F  are  the  chemical  symbols  for  the  carbon,  hydrogen,  and  fluorine 
atoms,  respectively  (figure  2).  The  carbon  atoms  are  connected  to  one  another  by 
single  bonds  to  form  the  polymer  backbone.  The  hydrogen  and  fluorine  atoms  are 
each  bonded  to  individual  carbon  atoms.  The  molecular  formula  of  PVDF  is  often 
written  as  (CHo-CFjjn,  where  n  typically  may  range  from  about  30,000  to  80,000.  The 
monomer  units  of  PVDF  can  be  linked  as  shown  in  figure  3  giving  rise  to  two  chain 
configurations.  The  CF2  portion  of  the  monomer  is  often  referred  to  as  the  "head" 
and  the  CH2  portion  the  "tail."  The  PVDF  chain  is  composed  largely  of  regular  head- 
to-tail  sequences.  Most  polymer  compositions  will  have  a  small  percentage  of  head- 
to-head  or  tail-to-tail  defects. 

P(VDF-TrFE)  is  synthesized  from  vinylidene  fluoride  and  1,1,2- 
trifluoroethylene.®  The  molecular  structure  of  the  copolymer  is  made  up  of  -CH2-CF2- 
and  -CHF-CF  2-  monomer  units  (figure  2)  connected  to  one  another  in  a  random 
fashion.  The  composition  of  a  P(VDF-TrFE)  copolymer  is  expressed  in  terms  of  the 
percent  distribution  of  moles  of  VDF  monomers,  e.  g.,  75  mol%  refers  to  a  copolymer 
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in  which  the  distribution  of  VDF  and  TrFE  monomer  units  is  3:1. 

In  the  melt  or  in  solution,  the  polymer  exists  as  numerous,  long,  spaghetti-like 
chains.  The  randomly  coiled  chains  are  constantly  moving.  These  motions  may 
consist  of  rotations  about  the  numerous  carbon-carbon  single  bonds  that  comprise  the 
polymer  backbone.^  During  cooling  of  the  melt  or  during  solvent  evaporation, 
dynamical  motions  about  the  carbon-carbon  single  bonds  lead  to  chain  folding.  The 
packing  of  the  folded  chains  results  in  the  formation  of  crystallites.  The  propensity 
of  these  polymer  chains  to  pack  into  discrete  crystal  structures  is  related  to  the  small 
van  der  Waals  (vdW)  radius  of  fluorine,  which  is  only  about  1.35  A.*  If  the  si2e  of 
fluorine  were  smaller,  let  us  say  about  the  si2e  of  the  hydrogen  atom  (1.2  A),  the 
fluoropolymers  would  behave  like  polyethylene  and  would  not  exhibit  as  many 
stable  polymorphic  structures  (the  structures  of  crystal  polymers  are  referred  to  as 
polymorphs).  The  polymer  would  then  tend  to  form  only  the  most  stable  poly¬ 
morph.  Were  the  van  der  Waals  radius  of  fluorine  larger,  the  rotational  barriers 
would  be  so  great  that  interconversions  between  polymorphs  would  be  restricted,  i.e., 
the  chain  would  stiffen.  This  is  not  the  case  with  fluoropolymers  where  the  barriers 
to  rotations  are  not  so  high  that  transitions  are  prevented  and  not  so  low  that  only 
the  most  stable  structure  is  always  favored.  In  other  words,  the  polymer  can  "lock" 
into  structures  not  necessarily  having  the  lowest  energies. 

The  crystal  structures  of  PVDF  and  P(VDF-TrFE)  can  be  described  as  stacks  of 
lamellae  immersed  in  a  sea  of  amorphous  regions  of  disordered  chain  conformations 
(figure  4).^'^°  The  lamellae  are  thin  platelet-like  structures  about  100  A  thick  and 
several  microns  broad. The  size  of  the  lamellae  is  sensitive  to  crystallization 
conditions.  A  microscopic  examination  of  the  interlamellar  regions  reveals 
nonordered,  isotropic  chains  connecting  the  crystallites.^'  The  chains  in  the  amor¬ 
phous  region  can  be  entangled  and  knotted.  Long  loops  can  exit  and  enter  the  same 
crystallite.  Some  loops  can  even  join  adjacent  lamellae.  Thus,  the  polymers  are 
complex,  two  phase,  semicrystalline  materials.  This  feature  of  PVDF  and  P(VDF- 
TrFE)  polymers  plays  a  profound  role  in  their  piezoelectric  performance. 

Although  several  polymorphs  for  PVDF  and  P(VDF-TrFE)  have  been  identi¬ 
fied,  we  will,  for  our  purposes,  consider  only  the  polymorphs  that  play  a  pivotal  role 
in  their  piezoelectricity,  the  a-  and  P-phases.  Segments  of  the  chain  conformations 
for  these  polymorphs  of  PVDF  are  illustrated  in  figure  5.  The  conformational 
structures  of  P(VDF-TrFE)  are  similar.  When  PVDF  crystallizes  from  the  melt  or 
solution  and  the  polymer  subsequently  anneals,  the  chains  pack  preferentially  into  a 
crystalline  polymorph  called  the  a-phase.’'^^  During  the  annealing  process,  the 
mobility  of  the  chains  guarantees  that  they  can  pack  into  crystal  cells.^^  In  the  a- 
phase,  the  sequential  arrangement  of  the  torsional  or  dihedral  angles  formed  between 
neighboring  C-F  and  C-H  bonds  as  one  looks  along  a  single  C-C  bond  is  in  the 
neighborhood  of  60°,  180°,  -60° ,  180°,  etc  (the  atoms  forming  a  dihedral  angle  are 
connected  to  one  another  sequentially  as  illustrated  in  figure  6).  For  shorthand,  we 
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write  this  as  tgtg'  where  t  refers  to  trans  or  180°,  and  g  or  gauche  is  a  60°  dihedral 
angle.  The  actual  torsional  angles  along  the  backbone  of  polymer  chain  may  deviate 
considerably  from  these  "ideal"  values.  However,  when  we  want  to  study  tfie  effects 
of  changes  in  dihedral  angles  on  chain  properties,  it  is  quite  useful  to  compare  the 
dihedral  angles  of  the  packed  chain  with  those  of  the  "ideal"  chain.  This  gives  us 
some  clue  concerning  the  magnitude  of  the  deformation  energies  required  to  distort 
the  chain  from  the  "ideal"  values.  This  loss  in  energy  stabilization  is  usually  compen¬ 
sated  for  by  a  great  gain  in  crystal  packing  energy.^^  This  is  an  important  point  to 
which  we  shall  return  below. 

One  of  the  salient  features  of  the  a-phase  crystal  is  that  it  is  nonpolar.  Two  a- 
phase  chains  pack  into  a  unit  cell  such  that  the  C-F  dipoles  normal  to  the  chain  axis 
cancel  out.  This  can  be  visualized  in  figure  7,  which  displays  a  projection  of  the  a- 
phase  unit  cell  looking  into  the  direction  of  the  chain  axis.  When  the  a-phase  is 
stretched  mechanically,  the  chains  undergo  torsional  motions  and  transform  into  a 
sequence  of  repeating  180°  torsional  angles.  This  is  the  P-phase,  which  is  character¬ 
ized  by  an  a\\-trans  conformation  (figures  5  and  7).  The  C-F  dipoles  of  the  3-phase 
chain  are  aligned  and  pointed  in  the  same  direction.  When  two  P-phase  chains  pack 
into  an  orthorhombic  unit  cell,  the  net  result  is  that  the  unit  cell  acquires  a  net  cell 
dipole  moment,  as  shown  in  figure  7.  When  poled  in  an  electric  field,  the  individual 
unit  cells  collectively  undergo  reorientation.  The  oriented  chains  are  now  polarized 
and  piezoelectric. 

P(VDF-TrFE)  is  known  to  crystallize  predominantly  into  a  P-phase  unit  cell 
possessing  monoclinic  symmetry.  The  net  polarization  of  a  unit  cell  containing  the 
CHF-CF,  is  less  than  that  of  CFj-CH,  since  the  former  has  a  C-F  dipole  on  CHF 
which  opposes  the  direction  of  the  dipole  on  CF,.  This  loss  in  polarization  might  be 
expected  to  result  in  the  loss  of  piezoelectricity,  but  this  is  not  so  as  shall  be  demon¬ 
strated  in  a  later  section.  We  shall  show  that  piezoelectricity  is  influenced  not  only 
by  the  polarization  of  the  crystalline  phase  but  also  by  the  mechanical  compliance  of 
both  the  crystalline  and  amorphous  phases.^^ 

The  fact  that  P(VDF-TrFE)  can  crystallize  into  the  polar  form  without  mechani¬ 
cal  stretching  is  one  of  the  great  advantages  of  the  copolymer.  However,  when  one 
considers  the  range  of  possible  P(VDF-TrFE)  compositions,  it  is  not  obvious  what 
ratios  of  VDF  and  TrFE  are  a  prerequisite  for  the  direct  crystallization  into  the  P- 
phase.  This  question  was  examined  from  the  perspective  of  molecular  modeling. 
Before  describing  this  study,  we  shall  provide  a  short  background  on  molecular 
modeling. 


4 


MOLECULAR  MODELING 


One  of  the  objectives  of  molecular  modeling  is  to  determine  the  optimum 
geometry  and  properties  of  a  molecule,  molecular  fragment,  or  crystal.  Many  useful 
properties  of  the  molecule  or  crystal  can  be  extracted  with  knowledge  of  the  opti¬ 
mized  geometry.  The  ability  to  predict  the  properties  of  solid-state  materials  using 
molecular  modeling  techniques  has  made  it  possible  for  material  scientists  to  tailor 
the  macroscopic  properties  of  these  materials  and  to  fine-tune  the  processing  condi¬ 
tions  used  to  engineer  these  materials  for  specific  applications.  Traditionally, 
molecular  modeling  was  used  by  chemists  to  study  tibe  properties  of  simple  organic 
and  inorganic  molecules  or  to  predict  the  course  of  chemical  reactions.  The  computa¬ 
tional  approaches  can  range  from  classical,  Newtonian  force  fields  that  treat  mole¬ 
cules  as  "ball  and  stick"  structures  to  quantum  mechanical  calculations  that  consider 
electronic  interactions.  The  approach  depends  on  the  nature  of  the  application  as 
well  as  the  size  of  the  molecular  structure.  One  can,  for  example,  solve  the 
Schrodinger  equation  for  the  hydrogen  molecule  with  a  simple  calculator.  But  for 
larger  molecular  systems  such  as  polymers,  adequate  computer  resources  are 
essential  to  calculate  certain  properties.  With  the  advent  of  powerful  computer  pro¬ 
grams  and  resources,  it  has  become  possible  to  extend  the  classical  and  quantum 
mechanical  algorithms  that  were  originally  intended  for  small  molecules  to  huge 
solid-state  systems  consisting  of  millions  of  atoms.  However,  even  with  large 
computers,  performing  calculations  on  solid-state  molecular  structures  is  still  a 
daunting  task.  A  far  better  approach  is  to  consider  fragments  of  the  solid-state.^^ 

The  task  can  be  made  even  more  manageable  if  the  solid-state  consists  of  periodic  or 
repeating  units  such  as  a  polymer.^®  By  focusing  on  the  properties  of  the  crystal  unit 
cell  and  the  dimensionality  of  the  crystal,  one  can  draw  conclusions  about  the 
behavior  and  characteristics  of  the  macroscopic  material.^® 

The  approach  in  this  study  was  to  perform  geometry  optimization  on  a  small 
segment  of  a  polymer  chain.“°  Geometry  optimization  and  assignment  of  charges 
were  carried  out  using  semiempirical  molecular  orbital  theory.  The  optimized 
structure  was  then  packed  into  a  crystal  cell.  Crystal  packing  calculations  were  then 
carried  out  on  the  relaxed  crystal  unit  cell  using  intermolecular  analytical  functions. 

In  calculating  the  energies  of  large,  single,  end-capped  polymer  chains,  we  used 
molecular  mechanics  for  economical  reasons  since  this  method  demands  fewer 
computer  resources  than  the  more  computationally  intensive  semiempirical  molecular 
orbital  theory.  Molecular  mechanics  was  also  used  to  calculate  the  net  dipole 
moments  of  crystal  unit  cells. 


MOLECULAR  MECHANICS 

In  the  molecular  mechanics  approach,  the  atoms  comprising  the  molecule  are 
visualized  as  hard  spheres  and  the  bonds  linking  the  atoms  as  springs.  The  energy 
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of  the  molecule  depends  on  the  interatomic  distances,  bond  angles,  conformations, 
and  relative  orientation  of  dipoles.  Each  interaction  can  be  described  by  an  analytical 
function  called  a  force  field.  From  the  geometry  of  the  equilibrium  structure,  it  is 
possible  to  deduce  certain  molecular  and  macroscopic  properties.  In  the  subsequent 
discussion,  it  will  be  useful  for  the  reader  to  refer  to  figure  6,  which  illustrates  the 
important  geometrical  features  of  a  molecule  including  the  bond  length,  bond  angle, 
torsional  or  dihedral  angle,  and  interatomic  distance. 

The  intramolecular  energy,  of  any  molecule  or  molecular  fragment  can  be 
written  as^^ 


^iiara  '^^sb  '^^vdw'^^tor 


(1) 

where  E^  is  the  bond  length  deformation  potential.  Eg  is  the  bond  angle  deformation 
potential,  E^j,  is  the  stretch-bend  potential,  E^.^w  is  the  van  der  Waals  interaction 
potential,  Ejo,  is  the  torsional  or  dihedral  potential,  and  E,,  is  the  dipole  interaction 
potential.  Each  term  in  equation  (1)  represents  a  component  of  the  total  force  field 
of  the  molecule  or  molecular  fragment.  The  force  field  due  to  the  sum  of  all  changes 
in  bond  lengths,  for  example,  can  be  described  by  Hooke's  law: 


(2) 


where  1  is  the  bond  length  between  two  atoms,  1^  is  the  equilibrium,  or  most  stable, 
bond  length,  and  is  the  force  constant  when  the  bond  length  is  deformed  or  dis¬ 
placed  from  the  lowest  energy  position.  The  bond-length  parameters  used  in  the 
molecular  mechanics  program  for  this  study  are  listed  in  table  1.^-'“ 

For  bond-angle  deformations,  the  potential  is 
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where  0  is  the  bond  angle,  0^  is  the  equilibrium  bond  angle,  and  kg  is  the  bond  angle 
deformation  force  constant.  The  relevant  parameters  are  shown  in  table  2.^'^ 
Stretch-angle  bending  cross  terms  refer  to  the  coupling  between  bond  stretching  and 
angle  bending.  These  motions  are  accommodated  by  the  expression  given  below: 


where  kji,  is  the  stretch-bending  force  constant.  The  stretch-bending  force  constants 
are  0.120  mdyn/  A  for  C-C-C,  C-C-F,  and  F-C-F  angles;  and  0.090  mdyn/  A  for  C-C-H 
and  H-C-F  angles.  No  force  constants  are  defined  for  H-C-H. 

The  analytical  function  for  torsional  or  periodic  rotations  about  single  bonds  is 


^tor=I^  [y  (1  -"y  (1  ‘^^^2(1))  3(1  +COS3<l))], 


(5) 


where  (j)  is  the  dihedral  angle  and  V^,  Vj,  V3  are  force  constants  for  individual 
torsional  barriers  (table  3).^ 

The  sum  of  non-bonded  van  der  Waals  interactions,  i.e.,  those  weak  forces 
arising  from  interactions  between  nonbonded  atoms  on  the  same  molecule,  can  be 
described  by 


^Hir=£€i,t2.9a:10^exp(-12.5p,y)-2.25p,y^, 


(6) 
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where  Sjj  is  defined  by 


and  Si  and  8j  are  referred  to  as  "hardness"  parameters.  The  following  expression 
defines 


(8) 


where  R^j  is  the  nonbonded  distance  between  the  interacting  atoms,  and  rjj  is  the  sum 
of  the  individual  van  der  Waals  radii  of  atoms  i  and  j  (r^j  =  rj  +  r^).  When  is  less 
than  3.311,  the  following  expression  is  used: 


^vdW  ^ 

The  nonbonded  parameters  used  in  equations  (6)  to  (7)  are  summarized  in  table  4. 
When  equation  (9)  is  used  for  a  C-H  bond,  s  and  r  are  0.046  and  3.340  A,  respective¬ 
ly- 

The  dipole-dipole  interactions  can  be  accommodated  by  the  following  equa¬ 
tion: 


_  [cosx-ScosajCosa; 
PiPj - T"^ - - 


where  and  Pj  are  the  dipole  moments  of  atoms  i  and  j;  is  the  dielectric  constant 
assumed  here  to  be  1.5;  %  is  the  angle  between  two  dipole  vectors;  and  oq  and  (x^  are 
the  angles  that  each  of  the  two  dipoles  make  with  the  vector  separating  them  (figure 
8). 

The  potential  energy  of  an  optimized  geometry  of  a  molecule  or  molecular 
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fragment  must  satisfy  the  equation 


(11) 


where  r;  represents  a  cartesian  coordinate.^  Typically,  the  coordinates  of  an  X-ray 
diffraction  structure  are  used  as  the  starting  or  initial  geometry.  The  first  derivative 
of  the  potential  energy  with  respect  to  the  atomic  coordinates  is  recorded  as  incre¬ 
ments  are  added  to  those  coordinates.  This  is  accomplished  by  changing  the  posi¬ 
tions  of  the  atoms  in  tiny  increments.  Minimization  proceeds  in  a  direction  of  the 
negative  gradient  of  the  potential  energy  until  equation  (11)  is  satisfied. 

The  success  of  a  molecular  mechanics  calculation  depends  on  the  selection  of 
realistic  force  constants.  These  are  generally  obtained  spectroscopically  or  by 
quantum  mechanical  calculations.  An  extensive  listing  of  published  force-field 
parameters  is  available."^ 

Geometry  optimizations  on  end-capped  polymers  were  carried  out  using  the 
force  fields  described  by  equations  (2)  to  (10).  Calculations  were  performed  using  the 
HyperChem  software  on  a  4D  Silicon  Graphics  workstation.  This  package  contains 
the  mm+  force  field,  which  is  derived  from  the  MM2  molecular  mechanics  software 
developed  by  Allinger."^  MM2  is  the  1977  version  of  molecular  mechanics,  which  is 
available  as  program  QCMP004  from  the  Quantum  Chemistry  Program  Exchange  at 
Indiana  University. 


SEMIEMPIRICAL  QUANTUM  MECHANICS 

Since  molecular  mechanics  does  not  explicitly  consider  effects  of  electrons,  the 
accuracy  of  this  methodology  depends  on  the  quality  of  the  force  fields  and  parame¬ 
ters  for  the  atoms  of  interest.  This  is  of  no  consequence  if  electronic  interactions  and 
bond-stretching  are  not  important.  To  obtain  "higher"  quality  results,  it  is  desirable 
to  consider  the  properties  of  the  electrons  and  to  solve  the  Schrodinger  equation  to 
obtain  these  properties.  This  procedure  can  be  done  with  semiempirical  quantum 
mechanics,  which  solves  an  approximate  form  of  the  Schrodinger  equation,  or  by  ab 
initio  quantum  mechanics,  which  does  not  use  any  experimental  parameters  to  solve 
the  Schrodinger  equation.  In  the  interest  of  economy  and  effciency,  we  chose  the 
semiempirical  approach  for  geometry  optimization  of  small  polymer  fragments.  The 
semiempirical  theory  was  also  used  to  calculate  the  tensile  elastic  moduli  when 
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stretching  polymer  chains.  This  operation  is  discussed  in  a  later  section. 

To  treat  electronic  effects  explicitly,  we  need  to  solve  the  Schrodinger  equation 
for  the  molecular  system  of  interest.  In  semiempirical  methods,  interactions  within 
the  nuclear  core  and  nonvalence  electrons  are  neglected.  Overlap  between  valence 
electrons  is  handled  by  empirical  functions. 

Semiempirical  methods  produce  approximate  solutions  to  the  Schrodinger 
equation: 


fr¥(r)=E^(r). 


(12) 


The  molecular  Hamiltonian,  H,  is  a  description  of  the  particles  of  a  system  and  'F(r) 
is  the  wave  function  of  the  molecular  system.  The  position  of  the  particles  of  the 
system  is  designated  by  r.  The  molecular  Hamiltonian  is  the  sum  of  kinetic  and 
potential  energy  terms: 


H=T+V, 


(13) 


where  T,  the  kinetic  energy,  is 


T= 


E— V, 


(14) 


where  mj,  is  the  mass  of  the  kth  particle.  The  differential  operator,  V,  is  given  by 


^  d  .  d  .  d  j 
V= — 1+ — 1+ — k, 
dx  dy  dz 


(15) 


where  x,  y,  and  z  represent  the  position  components  of  the  particle.  The  potential 
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energy,  V,  is 


i  I 
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u  ) 


(16) 


where  ARjj  is  the  distance  between  two  particles,  Z  is  the  atomic  number  of  the 
nucleus,  Ze  is  charge,  and  e  and  n  are  designations  for  electron  and  nucleus,  respec¬ 
tively.  The  lower  and  upper  case  subscripts  in  equation  (16)  also  refer  to  the 
elechons  and  nuclei,  respectively.  To  simplify  the  problem,  we  employ  the  Bom- 
Oppenheimer  approximation,  which  allows  us  to  treat  the  properties  of  nuclei  and 
electrons  separately.  The  Born-Oppenheimer  approximation  rests  on  the  fact  that 
nuclear  motion  is  much  slower  than  electronic  motion,  which  is  reasonable  given  that 
the  mass  of  an  electron  is  much  lighter  than  the  mass  of  the  nucleus.  Electronic 
motion,  which  can  now  be  treated  independently  without  considering  the  nuclear 
kinetic  energy,  is  described  by  the  electronic  Hamiltonian 
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When  the  Schrodinger  equation  is  solved  using  the  electron  Hamiltonian  from 
equation  (17),  the  energy  obtained  is  the  effective  nuclear  potential  function  describ¬ 
ing  the  potential  energy  surface  of  the  molecular  system.  It  is  important  to  note  that 
T  must  be  normalized;  that  is,  the  total  probability  of  finding  the  electron  in  space 
must  be 


f*“|cTpdv=l, 

*f  — oo 


(18) 


where  c  is  any  constant  and  v  is  "space  volume." 

There  are  a  number  of  simplying  assumptions  that  must  be  employed.  One  of 
these  is  that  the  molecular  orbital  can  be  expressed  as  a  linear  combination  of  the 
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individual  atomic  orbitals: 


(»’) 

where  (j)j^  is  an  atomic  orbital  and  C;^;  is  the  molecular  orbital  expansion  coefficient  for 
the  atomic  orbital.  By  convention,  we  designate  Roman  subscripts  for  molecular 
orbitals  and  Greek  subscripts  for  atomic  orbitals.  The  molecular  orbital  expansion 
coefficients  are  solved  by  using  the  variational  principle.  This  is  because  the  Hartree- 
Fock  method  is  based  on  the  relationship. 


£(S)>£(T), 


(20) 


i.e.,  the  minimum  energy,  E(H),  will  still  be  larger  than  the  exact  energy,  E('F). 

In  the  Hartree-Fock  method,  the  Self-Consistent  Field  (SCF)  approach  is 
employed  where  the  position  of  each  electron  is  guessed  as  well  as  its  average 
potential  in  the  field  of  other  electrons.  The  Schrodinger  equation  is  solved  to  obtain 
a  new  "guess"  of  the  positions.  This  procedure  is  repeated  until  the  wavefunction  of 
each  electron  is  "consistent"  with  the  fields  of  all  electrons.  When  solving  for  the 
coefficients  in  equation  (19),  the  energy  of  the  wavefunction  must  be  minimized.  To 
solve  for  the  molecular  orbitals,  we  must  determine  the  coefficients  of  the  atomic 
orbitals.  We  do  this  step  using  the  Roothaan-Hall  secular  equation: 

1 

where  e-  is  the  spin-free  one-electron  orbital  energy  of  an  orbital  (ji  =  1,  2,  •••,  N);  S  is 
the  matrix  representing  overlap  between  orbitals,  i.e., 

5,v  =  ^9,(l)|<Pv(l)^  (22) 


and  F  is  the  Fock  matrix. 


N  N 

k  a 

where  is  a  matrix  of  the  energy  of  a  single  electron  in  the  field  of  bare  nuclei. 
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The  term  ^jiv  |  >  is  a  two-center  two-electron  integral  between  a  pair  of  atoms;  thus. 


^|iv|Ao>=j'(p^(l)(p^(l) 


'  1^ 


V  12/ 


‘Pi(2)(Pa(2)i/t, 


(24) 


where  is  the  density  matrix,  which  can  be  written  as 


occtq>ied 

i=l 

Solving  equation  (21)  is  done  by  matrix  diagonalization.  This  leads  to  the  coefficients 
used  to  describe  the  atomic  orbitals  in  equation  (19).  The  total  electronic  energy  in 
equation  (21)  is,  therefore,  defined  as 


E=^P{H+F).  (26) 

The  details  of  matrix  diagonalization  and  the  functions  used  to  calculate  molecular 
properties  such  as  atomic  charge  are  outside  the  scope  of  this  report.  The  interested 
reader  may  want  to  consult  one  of  the  excellent  reviews  or  books  referenced  in  this 
section. 

MOP  AC  is  currently  the  most  popular  general-purpose  semiempirical  molecu¬ 
lar  orbital  program.^^  This  program  represents  the  collective  effort  of  numerous 
groups  who  over  the  years  have  implemented  a  variety  of  semiempirical  methods.  In 
MOP  AC,  a  starting  geometry  must  be  specified.  As  in  the  case  of  classical  methods, 
geometry  optimization  is  carried  out  with  derivatives  with  respect  to  Ihe  coordinates, 
hi  each  iteration,  the  heat  of  formation  is  calculated.  When  there  is  no  further  change 
in  energy,  geometry  optimization  stops.  In  this  work,  MOP  AC  was  used  to  minimize 
the  polymer  segments  which  were  packed  into  crystal  cells.^^  This  was  done  to 
accommodate  the  effects  of  the  atomic  charges  on  the  crystal  packing  energies. 


CRYSTAL  PACKING  CALCULATIONS 

The  structure  obtained  after  minimization  of  a  chain  segment  is  used  to  calcu¬ 
late  the  total  potential  energies  and  the  dipole  moments  of  individual  unit  cells.  To 
perform  this  operation,  it  is  first  necessary  to  "pack"  the  individual  chain  segments 
into  appropriate  unit  cells.  "Packing"  refers  to  the  alignment  and  arrangement  of  the 
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chain  segments  into  a  crystal  cell.  The  lattice  dimensions  and  angles  are  altered  to 
minimize  repulsive  electrostatic  and  vdW  interactions  between  atoms  on  neighboring 
chains  and  to  maximize  attractive  electrostatic  and  vdW  interactions.  Intermolecular 
electrostatic  interactions  are  accommodated  by  the  sum  of  individual  interactions 
between  charged  atoms  on  neighboring  chains:^ 


where  is  the  coulombic  energy,  and  q  are  the  charges  on  atoms  i  and  j.  Atoms 
possessing  the  same  charge  will  repel  one  another,  whereas  those  with  opposite 
charges  attract  Intermolecular  van  der  Waals  interactions  are  calculated  using  the 
potential  shown  below:^ 


(28) 

One  usually  starts  with  the  structure  of  a  crystal  whose  dimensions  have  been 
determined  by  x-ray  crystallography.  Once  the  chain  segments  are  packed  within  the 
cell,  the  interactions  between  them  are  minimized  by  changing  the  lattice  dimensions 
and  angles  in  tiny  increments.  Eventually,  a  relaxed  crystal  structure  and  its  intermo¬ 
lecular  energy  are  obtained.  The  geometry  of  this  structure  should  not  deviate 
significantly  from  the  experimental  x-ray  structure.  The  crystal  packing  calculations 
described  in  this  report  were  performed  using  the  CERIUS  Version  3.1  software 
provided  by  Molecular  Simulations. 


RESULTS  AND  DISCUSSION 

Before  designing  a  transducer,  it  is  desirable  that  the  properties  of  the  active  or 
passive  piezoelectric  material  be  fully  characterized.  In  the  case  of  P(VDF-TrFE), 
which  is  the  material  of  choice  for  an  application  requiring  thick  slabs  or  unusually 
shaped  designs,  one  needs  first  to  select  a  suitable  ratio  of  VDF  and  TrFE  that 
crystallizes  directly  into  the  desired  P-phase  and  that  satisfies  performance  specifica¬ 
tions  of  the  transducer  application.  Thus,  one  of  our  primary  objectives  was  to  deter¬ 
mine  the  effect  of  fluorine  content  on  the  relative  stabilities  of  the  a-  and  P-phases  of 
PVDF  and  P(VDF-TrFE).  We  will  discuss  in  the  following  sections  how  molecular 
modeling  was  used  to  study  (1)  the  effect  of  polymer  composition  of  the  crystalliza¬ 
tion  behavior  of  PVDF  and  P(VDF-TrFE)  and  (2)  the  determination  of  the  electrical 
and  mechanical  properties  of  these  materials. 
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THE  CRYSTALLIZATION  BEHAVIOR  OF  PVDF  AND  P(VDF-TrFE) 


The  crystal  structures  favored  during  crystallization  of  polymer  melts  or  solu¬ 
tions  are  determined  by  the  combined  intermolecular  and  intramolecular  energies. 
The  intermolecular  energies  refer  to  the  van  der  Waals  and  Coulombic  interactions 
between  neighboring  chains  packed  in  the  unit  cells.  The  intramolecular  energies  are 
determined  by  the  chain  conformation,  i.e.,  the  number  and  position  of  the  trans  and 
gauche  linkages.  The  procedure  used  to  study  the  energetics  of  crystallization  is 
outlined  in  figure  9. 

The  initial  structures  used  to  build  up  larger  chain  segments  were  dimeric 
segments.  The  atomic  positions  of  the  dimers  were  determined  by  Hasegawa  et  al., 
who  employed  x-ray  crystallography.^  These  workers  proposed  that  the  PVDF  a- 
phase  is  a  distorted  tgtg'  structure,  with  a  carbon  dihedral  repeating  pattern  of  179° 
and  45°.^  According  to  the  Hasegawa  study,  the  P-phase  conformation  is  a  slightly 
distorted  aW-trans  alternating  planar  zig-zag  with  repeating  172°  and  -172°  angles. 

P(VDF-TrFE)  polymer  segments  were  constructed  from  the  PVDF  segment  by 
replacing  the  appropriate  number  of  hydrogen  atoms  with  fluorine  atoms.  The  bond 
lengths  and  bond  angles  of  the  Hasegawa  a-  and  3-phase  structures  are  shown  in 
table  5.  Dimeric  structures  of  PVDF  and  P(VDF-TrFE)  were  first  minimized  using 
semiempirical  quantum  mechanics.  Atomic  charges  were  calculated  and  assigned  to 
each  atom  (figure  10).  These  segments  were  linked  to  together  to  form  chains 
consisting  of  sixteen  end-capped  carbon  atom  structures.  An  extended  structure  of  a 
PVDF  chain  is  depicted  in  figure  11.  Only  one  chain  of  PVDF  was  considered,  since 
only  one  stereoregular  arrangement  of  this  polymer  is  possible.  Although  the 
structure  of  the  carbon  backbone  of  P(VDF-TrFE)  is  similar  to  that  shown  in  figure  11 
(with  the  exception  that  some  hydrogen  atoms  are  replaced  by  fluorine  atoms), 
several  arrangements  of  P(VDF-TrFE)  are  possible,  however,  since  the  monomers  of 
these  polymers  can  adopt  isotactic,  syndiotactic,  and  atactic  configurations  (figure  12). 
By  configuration,  we  refer  to  the  arrangement  of  fluorine  and  hydrogen  atoms  on 
each  carbon  center.  In  an  isotactic  polymer,  all  of  the  monomers  have  the  same 
configuration.  In  an  syndiotactic  polymer,  the  configuration  alternates  from  one 
monomer  to  the  next.  Atactic  structures  are  random.  In  this  study,  only  isotactic 
and  syndiotactic  structures  were  considered.  Twelve  structures  of  75  mol%  P(VDF- 
TrFE),  where  two  TrFE  monomers  are  separated  by  three  VDF  monomer  units,  were 
selected.  The  configurations  of  the  two  -CFH-  centers  represent  isotactic  and 
syndiotactic  structures.  Each  structure  of  the  50  mol%  P(VDF-TrFE)  composition  was 
also  selected  to  represent  isotactic  and  syndiotactic  configurations.  The  sequences  of 
all  polymer  chains  consisted  of  head-to-tail  linkages.  No  head-to-head  or  tail-to-tail 
sequences  were  considered. 

The  bond  lengths,  bond  angles,  and  dihedral  angles  of  the  single  carbon  end- 
capped  chains  were  optimized  using  molecular  mechanics.  The  minimizations  were 
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performed  by  allowing  all  of  the  atoms  in  the  chain  to  relax.  The  resulting  conforma¬ 
tions  were  helical  structures.  The  energies  of  these  structures,  which  are  tabulated  in 
table  6,  increase  with  fluorine  content.  The  Hasegawa  and  "ideal"  structures  were 
then  minimized  via  relaxation  of  the  hydrogen  and  fluorine  atoms  while  the  positions 
of  the  carbon  atoms  in  the  polymer  backbone  chain  were  fixed.  This  approach 
allowed  us  to  preserve  the  original  geometry  of  the  experimental  chain  conformation. 
When  we  compare  the  energies  of  each  structure,  we  note  that  the  Hasegawa  fgtg' 
structure  is  more  stable  than  the  Hasegawa  and  "ideal"  all-trans  conformations.  But 
the  "ideal"  all-trans  structure  is  favored  for  the  copolymers.  A  plot  of  the  data  (figure 
13)  shows  that  the  crossover  point,  i.e.,  the  composition  where  the  "ideal"  all-frans 
becomes  more  stable  than  the  tgtg'  is  about  90  mol%  VDF.  These  calculations  agree 
with  the  well-documented  experimental  observation  that  PVDF  is  favored  to  crystal¬ 
lize  into  the  tgtg'  conformation,  whereas  the  copolymer  crystallizes  into  the  all-frans 
structure.^^  The  plots  in  figure  13  underscore  the  small  differences  in  energy 
between  the  conformations.  During  annealing,  however,  formation  of  the  most 
stable  conformations  will  be  favored. 

To  build  crystal  unit  cells  of  PVDF  and  P(VDF-TrFE),  chain  segments  were 
packed  into  a  crystal  unit  cell.  The  direction  of  the  polymer  chain  was  aligned  along 
the  c-axis,  and  the  structure  was  subsequently  centered.  There  were  no  symmetry 
restrictions  imposed  on  each  unit  cell.  Starting  with  the  initial  structures  based  on  x- 
ray  data,  the  dimensions  of  the  unit  cells  were  minimized  iteratively  by  calculating 
the  energies  of  minor  deformations  in  the  lattice  structures  as  well  as  rotations  of  5\e 
polymer  chains  about  the  c-axis  (figure  14).  The  relaxed  cell  parameters  of  the  PVDF 
and  P(VDF-TrFE)  unit  cells  are  summarized  in  table  7.  In  agreement  with  x-ray 
studies,  the  a  and  b  crystal  dimensions  increase  with  increasing  fluorine  substitution. 
This  expansion  of  the  cell  lattice  dimensions  can  be  attributed  to  the  increasing  van 
der  Waals  energies  as  more  fluorine  atoms  are  added  to  the  crystal  structure. 
Examining  the  crystal  energies  (table  8),  we  note  that  the  Hasegawa  tgtg'  crystal 
structure  is  favored  over  the  "ideal"  structure  for  all  compositions.  This  calculation 
supports  the  Hasegawa  structure,^^  and  suggests  that  during  crystallization  from  the 
melt  or  solution,  the  t  and  g  dihedral  bonds  suffer  distortions  as  the  chains  within  the 
crystal  lattice  rotate.  The  energy  required  to  distort  the  dihedral  angles  is  probably 
more  than  offset  by  the  energy  gained  during  crystal  packing. 

The  intramolecular  energies  calculated  on  the  end-capped  polymer  chain 
segments  and  the  resulting  crystal  packing  energies  are  not  additive  in  this  study 
since  these  models  do  not  have  the  same  number  of  atoms.  However,  we  can 
compare  these  energies  and  establish  several  important  trends.  The  crystal  packing 
energies  of  PVDF  and  P(VDF-TrFE)  favor  the  all-fra?is  conformation.  For  the  copoly¬ 
mers,  this  is  in  line  with  the  propensity  of  these  compositions  to  crystallize  into  the 
all-trans  3-phase.  But  the  crystal  packing  calculations  contradict  the  crystallization 
behavior  of  PVDF,  which  is  known  to  crystallize  directly  into  the  tgtg'  a-  phase. 

One  explanation  for  this  discrepancy  is  that  the  intramolecular  interactions  play  a 


16 


vital  role  in  controlling  the  kinetics  of  crystallization.  During  crystallization,  although 
the  energy  of  the  a\l-trans  PVDF  unit  cell  is  lower  than  the  energy  of  the  fgfg'  cell,  a 
large  energy  barrier  may  prevent  crystallization  into  the  P-phase.  This  barrier  is 
somewhat  relaxed  for  the  copolymer  compositions  so  that  crystallization  into  the  3- 
phase  is  favored  both  kinetically  and  thermodynamically.  Recent  calculations  using 
molecular  dynamics  have  been  undertaken  in  this  laboratory  and  support  this  idea. 
These  results  will  be  reported  in  the  near  future.^^ 

The  most  important  conclusion  of  this  portion  of  the  study  is  that  these 
calculations  demonstrate  that  crystallization  into  the  desired  3-phase  is  favored  for 
P(VDF-TrFE)  compositions  where  VDF  is  less  than  90  mol%.  This  knowledge 
allowed  us  to  set  an  upper  limit  on  the  percent  composition  of  VDF  to  achieve  the  3- 
phase  without  mechanical  stretching.  Ideally,  the  VDF  concentration  should  be  as 
high  as  possible  to  also  maximize  the  Curie  temperature,  which  is  the  crystal-crystal 
transition  between  the  lower  temperature  ferroelectric  phase  and  the  higher  tempera¬ 
ture  paraelectric  phase.  Above  the  Curie  temperature,  the  polymer  looses  its 
polarization.  Since  the  Curie  temperature  increases  from  about  125°C  at  VDF  w  60 
mol%  to  over  180°C  for  PVDF,^'^*'^^  to  satisfy  the  contradictory  requirements  of 
energetics  and  temperature  stability,  the  most  suitable  copolymer  compositions 
should  be  those  with  VDF  content  between  60  and  90  mol%. 


CALCULATION  OF  THE  ELASTIC  MODULI  AND  PIEZOELECTRIC 
CONSTANTS 

When  a  piezoelectric  polymer  unit  cell  is  stressed,  the  resulting  strains  along 
the  c-axis  (the  chain  direction),  a-axis,  and  b-axis  (poling  direction)  can  be  written  as 


and 


5c  -Si  1  Tj  +S12T2  +Si2,Ty 

(29) 

~^21  ^1  ■'■^22^2  ■'■‘^23^3’ 

(30) 

“^31  ■‘■'^32^2  ■‘■■^33^3’ 

(31) 

where  S^,  S^,  and  are,  respectively,  the  strains  along  the  c,  a,  and  b  directions,  Sjj'^  is 
the  real  elastic  compliance  (i,j  =  1,  2,  and  3),  and  T^,  Tj,  and  T3  are  the  applied 
stresses  in  the  c,  a,  and  b  directions,  respectively  (figure  15).  When  a  stress,  T3,  is 
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applied  to  an  undamped  crystal  cell  along  the  direction  of  the  poling  field,  then 


(S33 +523 +Si3)7’3  ^33^3* 


(32) 


where  833*^'  is  the  compliance  of  an  undamped  crystal  cell.  We  will  assume  for  now 
that  the  c-axis  is  clamped,  that  is  when  the  b-axis  is  stretched,  the  a-axis  deforms  but 
the  c-axis  does  not  change  its  dimension.  Thus,  from  equation  (29),  the  stress  of  a 
crystal  clamped  in  the  chain  or  c-axis  direction  (S^  =  0;  T,  =  0)  can  be  written  as 


(33) 


Note  that  we  have  assumed  that  8^3*^  =  835*^.  So  if  the  c-axis  is  clamped  and  =  0, 
then 


S/j  +S^  -(521  +5'3i)7’j  +(523 


(34) 


With  the  further  assumption  that  83^*^  «  s,!*",  we  obtain  from  equations  (33)  and  (34) 


2(53^)^ 


[533  +523  - : — J  ^3  --S33 


T 

I3. 


Ml 


(35) 


The  effective  compliance  of  the  clamped  cell,  S33'",  is  therefore 


(36) 


The  procedure  to  calculate  the  elasHc  compliances  is  based  on  a  curve-fitting 
method  to  calculate  the  tensile  compliance  of  long-chain  polymers  with  molecular 
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modeling  techniques.^®  As  adapted  to  the  present  work,  the  procedure  consisted  of 
the  application  of  tiny  homogeneous  stresses  to  the  plane  of  the  polymer  unit  cell 
normal  to  the  poling  or  "3"  direction.  In  the  case  of  a  PVDF  crystal,  since  the  unit 
cell  is  orthorhombic,  only  the  a-dimension  was  permitted  to  undergo  relaxation.  The 
c-dimension  was  clamped.  The  intermolecular  energy  consisting  of  van  der  Waals 
and  electrostatic  interactions  was  calculated  for  each  deformed  structure.  The 
energies  and  deformations  were  then  subjected  to  curve-fitting  by  linear  regression  to 
the  cubic  form  of 


(37) 


where  K  is  the  force  constant  of  the  deformation  Ab,  +  e  the  vdW 

and  electrostatic  energies,  and  A,  B,  and  C  are  constants.  The  force  constant,  K,  was 
calculated  from  the  slopes  of  equation  (37).  The  correlation  coefficient  of  each  curve 
was  at  least  0.9998.  As  exemplified  by  a  plot  of  K  vs.  elongation  of  a  PVDF  unit  cell 
(figure  16),  the  force  constant  typically  displayed  wide  excursions  at  deformations 
less  than  0.05%.  Accordingly,  the  force  constants  were  determined  within  the  ranges 
of  deformations  represented  by  the  flat  regions  of  the  curves.  Thus,  the  elastic 
compliance  associated  with  changes  in  the  b-dimension  of  the  crystal  unit  cell,  is 


c''_  1  ac 

‘Tb' 


(38) 


where  ac  is  the  area  of  the  plane  perpendicular  to  the  direction  of  the  deformation. 

The  elastic  compliance,  s^*^,  was  calculated  by  simulating  a  stretching  of  a 
single  polymer  chain  along  the  c-axis.^^  A  semiempirical  quantum  mechanical 
approach  was  chosen  for  chain  minimization  since  a  previous  semiempirical  study 
demonstrated  good  agreement  between  calculated  and  experiment  tensile  moduli  of 
several  polymer  chains.^®'^"  In  this  case,  the  heat  of  formation  of  each  structure  was 
calculated.  Each  of  these  chain  structures  was  then  packed  into  a  crystal  cell  whose 
a-  and  b-dimensions  were  subsequently  minimized.  The  minimized  energies  and  b- 
dimensions  were  then  used  to  calculate  s^g*^. 

The  elastic  compliances  are  tabulated  in  table  9.  The  values  listed  here  show 
that  the  difference  between  Sgg'^'  and  Sgg'"  of  a  unit  cell  of  PVDF  and  the  50  mol% 
copolymer  is  less  than  1%,  i.  e.,  s^'  «  Sgg*^".  This  can  also  be  assumed  to  be  the  case 
for  the  75  mol%  P(VDF-TrFE). 

We  now  turn  to  the  calculation  of  changes  in  polarization  when  a  crystal  cell  is 
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deformed.  Stressing  a  polymer  unit  cell  in  the  direction  of  the  poling  field  induces  a 
change  in  its  volume  and  saturation  polarization:^ 


p/_p ' 
*  S  ^  S 


(39) 


where  dgj"  is  the  piezoelectric  constant,  and  P/  and  P/  are  the  initial  and  final 
saturation  polarizations.  We  define  the  saturation  polarization  of  a  polymer  unit  cell 

^g43, 44,45 


(40) 


where  jii.  is  the  dipole  moment  of  an  individual  C-F  bond  and  v,  is  the  volume  of  the 
unit  cell.  The  procedure  to  calculate  the  total  dipole  moments  of  PVDF  and  P(VDF- 
TrFE)  crystal  cells  is  based  on  the  rigid  dipole  model.^^  In  this  approach,  the  dipole 
moments  of  the  individual  C-F  within  a  packed  unit  cell  were  summed  to  obtain  the 
total  dipole  moment  of  the  cell.  The  changes  in  cell  volume  and  polarization,  and  the 
piezoelectric  constants  were  calculated  by  simulating  stresses  along  the  b-axis  of  the 
minimized  crystal  structures  (table  9).  The  calculated  Sjs^'  elastic  compliance  and  djj" 
piezoelectric  constants  of  the  ciy^stals  are  also  tabulated  in  table  9.  These  data 

that  the  crystals  of  PVDF  and  P(VDF-TrFE)  are  "hard"  materials  and  that 
■fluorine  substitution  decreases  the  compliance  along  the  b-axis.  A  possible  explana¬ 
tion  for  the  low  compliance  may  be  found  in  the  interactions  taking  place  during 
deformation  of  the  lattice  dimensions.  Changes  in  the  lattice  geometry  of  the  cell  are 
accompanied  by  unfavorable  van  der  Waals  interactions  between  fluorine  atoms  on 
adjacent  chains. 

The  values  here,  which  are  smaller  than  those  obtained  by  Goddard  and 
Karasawa'*^  and  Tashiro  et  al./^  suggest  that  polymer  crystals  by  themselves  are  low 
compliance  materials  and  therefore  possess  low  levels  of  piezoelectricity.  Flowever, 
the  results  of  the  calculations  in  this  work  are  not  inconsistent  with  the  well-known 
piezoelectric  properties  of  semicrystalline  PVDF  and  P(VDF-TrFE),  if,  as  will  be 
demonstrated  below,  the  origin  of  piezoelectricity  in  the  semicrystalline  polymers  lies 
in  the  properties  of  both  the  crystalline  and  amorphous  phases.  In  this  model. 
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(41) 


tiie  Young's  modulus  of  a  semicrystalline  polymer  can  be  written  as^ 

where  ¥33®'^  and  Y33®  are  the  Young's  moduli  of  the  semicrystalline  polymer  and  amor¬ 
phous  phase,  respectively,  and  X<.  is  the  degree  of  crystallinity.  The  Young's  modulus 
of  a  semicrystalline  polymer  of  known  crystallinity,  Y  33“,  can  be  obtained  from  the 
literature.  On  substitution  of  the  literature  value  into  equation  (41),  S33‘'  can  be 
calculated  (Y33  =  I/S33).  It  is  likely  that  S33®  is  independent  of  crystallinity,  so  that 
equation  (41)  can  also  be  used  to  estimate  S33“  of  piezoelectric  polymers  possessing 
varying  crystallinities. 

The  piezoelectric  constant  of  a  semicrystalline  polymer  is  calculated  from^ 


(«) 

where  d33‘^  is  the  piezoelectric  response  due  to  dimensional  changes,  and  0  is  the 
average  orientation  of  dipoles  along  the  poling  direction. 

We  have  calculated  the  S33  elastic  compliances  and  d33  constants  for  hypotheti¬ 
cal  semicrystalline  PVDF  and  75  mol%  P(VDF-TrFE)  compositions  using  equations 
(41)  and  (42).  These  are  the  predicted  values  for  semicrystalline  polymers.  Consider 
first  the  calculation  for  PVDF.  An  experimental  compliance  of  PVDF  with  =  0.5 
was  taken  to  be  4.72  x  10'^°  m“/N.^^  Elastic  compliance  values  were  calculated  for 
PVDF  possessing  different  degrees  of  crystallinity.  From  the  elastic  compliance  and 
polarization,  we  calculated  the  piezoelectric  constant.  As  shown  in  table  10,  the 
experimental  d33  values  agree  well  with  the  calculated  d33‘^  values  for  those  composi¬ 
tions  where  X^  =  0.5.  For  75  mol%  P(VDF-TrFE),  S33®  was  calculated  using  the  experi¬ 
mental  elastic  compliance,  3.0  x  10‘^°  m^N  for  =  0.75}^  We  took  into  account  the 
possibility  of  local-field  effects  for  75  mol%  P(VDF-TrFE).  It  is  known  that  local-field 
effects  may  enhance  the  saturation  polarization  of  the  75  mol%  P(VDF-TrFE)  unit 
cell.^  Including  this  local-field  correction  improves  the  agreement  between  the 
experimental  and  calculated  djj'*  constants  for  the  75  mol%  P(VDF-TrFE)  copolymer 
(local-field  corrections  for  PVDF  can  be  neglected^'^^).  Again,  as  in  the  case  of  PVDF, 
the  agreement  between  the  experimental  and  calculated  piezelectric  constant  is  good. 

There  are  several  conclusions  to  be  drawn  from  this  treatment.  First,  the 
maximum  dgj*^  occurs  at  compositions  that  are  about  50%  crystalline,  and  d33‘^  is 
predicted  to  decrease  at  higher  crystallinities.  This  is  not  surprising  in  light  of  our 
calculation  that  predicts  "hard"  crystallites.  Another  conclusion  is  that  a  significant 
decrease  in  piezoelectricity  should  occur  at  low  crystallinites  due  to  low  crystal 
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polarizations.  Both  of  these  conclusions  allow  us  to  portray  the  microstructure  of 
semicrystalline  polymers  as  a  collection  of  oriented  C-F  dipoles  embedded  in  a  "soft" 
amorphous  phase.^"  The  piezoelectric  response  is  therefore  a  result  of  small  motions 
of  these  dipoles.  The  key  feature  of  this  model  is  that  the  polymer  can  be  treated  as 
a  composite  of  tiny  but  polarized  crystallites  immersed  in  a  nonpolarized,  relatively 
soft  amorphous  polymer. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  original  objective  of  this  work  was  to  demonstrate  that  molecular  model¬ 
ing  is  a  practical  methodology  to  aid  in  the  proper  selection,  preparation,  and 
processing  of  piezoelectric  fluoropolymers  for  transducer  applications.  The  impetus 
for  such  an  undertaking  was  to  broaden  our  practical  knowledge  of  these  complex 
materials.  It  was  hoped  that  by  modeling  the  atomic  and  molecular  properties  of 
PVDF  and  P(VDF-TrFE)  in  an  effort  to  identify  the  most  suitable  compositions  for 
hydrophone  and  projector  applications  that  sufficient  "in-house"  knowledge  would  be 
acquired  to  maximize  the  practical  application  of  these  materials.  Such  an  effort 
would  naturally  lead  to  a  procedure  where  one  would  select  the  composition  and 
subsequent  processing  conditions  from  a  knowledge  of  the  final  transducer  properties 
and  performance.  This  procedure  might  be  called  "molecular  tailoring"  or  design  of 
transducer  materials,  and  could  be  conceived  as  an  algorithm  where  the  transducer 
designer  might  enter  the  desired  parameters  of  a  transducer  application  such  as  a 
figure  of  merit  (FOM)  and  obtain  the  materials  properties  such  as  degree  of 
crystallinity,  polarization,  and  elastic  complicance  diat  would  result  in  the  desired 
FOM.  That  we  are  at  the  threshold  of  such  a  procedure  is  evidenced  by  our  success 
in  confirming  certain  well-known  properties  of  PVDF  and  the  copolymers  through 
computational  chemistry  techniques  and  by  the  success  we  have  had  in  building  a 
multilayered  copolymer  stack  for  an  active  projector.^  In  the  meantime,  this  project 
has  lead  to  several  interesting  and  practical  conclusions  concerning  PVDF  and 
P(VDF-TrFE). 

1.  We  have  identified  P(VDF-TrFE)  compositions  that  should  be  suitable  for  most 
piezoelectric  applications  requiring  unusual  shapes  or  thick  slabs  of  copolymer. 

These  may  be  applications  where  mechanical  stretching  (as  in  the  case  of  PVDF) 
would  be  difficult  to  perform.  The  calculations  in  this  work,  which  were  carried  out 
on  copolymer  model  molecules,  suggest  that  compositions  containing  between  60  and 
90  mol%  VDF  crystallize  directly  into  the  |3-phase,  have  sufficient  polarization,  and 
possess  Curie  transitions  sufficiently  high  to  rule  out  the  possibility  of  depoling 
during  long  periods  of  use  under  moderate  temperature  conditions.  It  is  recom¬ 
mended  that  attention  be  given  to  those  copolymers  having  these  levels  of  VDF. 

2.  The  results  of  the  elastic  compliance  study,  although  still  tentative,  suggest  that 
piezoelectric  crystals  are  by  themselves  low  compliant  materials  and  that  piezoelec- 
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tricity  arises  from  the  collective  motion  of  dipoles  embedded  in  the  soft  amorphous 
phase  of  the  polymer.  If  this  is  true,  then  the  optimum  crystallinity  would  be  at 
about  Xc  =  0.5.  This  degree  of  crystallinity  can  be  achieved  by  carefully  monitoring 
the  annealing  conditions  and  by  performing  differential  scanning  calorimetry.  At  this 
stage,  it  is  recommended  that  more  attention  be  given  to  the  effects  of  crystallinity  on 
piezoelectric  performance  and  also  on  long  term  aging.  It  is  also  suggested  that 
further  investigation  be  given  to  checking  and  testing  the  hypothesis  suggested  by 
the  results  of  this  study  that  piezoelectric  fluoropolymeric  crystals  are  indeed  "hard." 

3.  A  final  recommendation  is  suggested  by  the  focus,  in  this  study,  on  the  molecular 
properties  of  fluoropolymers.  Although  attention  was  given  to  PVDF  and  P(VDF- 
TrFE),  there  are  other  molecular  combinations  to  consider.  With  only  a  few  excep¬ 
tions,  little  work  has  been  done  on  other  fluoropolymeric  systems,  including,  for 
example,  copolymers  of  vinyl  fluoride  and  vinylidene  fluoride,  tetrafluoroethylene 
and  vinylidene  fluoride,  to  name  only  a  few.  It  would  be  a  fruitful  enterprise  to 
study  the  intramolecular  and  crystal  packing  energies  of  these  copolymers  as  well  as 
their  polarizatioris  and  elastic  compliances,  using  the  techniques  described  in  this 
report,  and  to  compare  these  values  with  the  data  for  PVDF  and  P(VDF-TrFE).  Such 
an  enterprise  might  reveal  superior  performance  properties  and  lead  to  the  develop¬ 
ment  of  new  and  exotic  transducer  materials. 

In  conclusion,  this  project  demonstrates  the  validity  and  practicality  of  using  a 
"molecular"  approach  for  the  study  of  material  properties.  As  molecular  modeling 
matures,  computational  chemistry  will  emerge  as  a  powerful  tool  to  sort  out  the 
complexities  of  candidate  materials  for  practical  applications  such  as  described  in  this 
report 
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HOMOPOLYMER; 


-  A- A- A-A- A- A- 


COPOLYMER: 


-A-B-A-B-A-B-  (ALTERNATING) 
-A-A-B-A-B-B-A-B  -  (RANDOM) 
A  ,  B  =  MONOMERS 


Figure  1.  Representation  of  the  Chain  Structures  of  Polymers 


PVDF: 


H  F  H  F 

I  I  I  I 

A/V- 0  -  c  -  c  -  c  VV\ 

I  I  i  I 

H  F  H  F 


P(VDF  -TrFE): 


H  F  H  F 

I  I  I  I 

AA^  c  -  c  -  c  -  c  VsA 
.1111 
F  F  H  F 


Figure  2.  Extended  Molecular  Structures  of  PVDF  and  P(VDF-TrFE)  Chain  Segments 
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Figure  3.  Extended  Molecular  Structures  Showing  Chain  Connectivity 
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CRYSTAL 

REGION 


FOLDED 

LAMELLAE 


FROM:  STEVENS,  "POLYMER  CHEMISTRY:  AN  INTRODUCTION", 

2ND  ED.,  OXFORD  UNIVERSITY  PRESS:  OXFORD,  1990,  P.92 


Figure  4.  Stacks  of  Lamellar  Polymers 
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HYDROGEN 

CARBON 

FLUORINE 


a  PHASE 


P  PHASE 


HYDROGEN 

1  FLUORINE 

Figure  7.  Crystal  Structures  of  the  a-  and  P-Phases 
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Figure  9.  Flow  Chart  Used  to  Study  the  Energies  of  Crystals  and  Chains 
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PVDFTT 


PVDF  TGTG 


P(VDF-TrFE)  TT 


Figure  10.  Atomic  Charges  Calculated  by  Semiempirical  Quantum  Mechanics 


Figure  11.  Extended  Structure  of  PVDF  End-Capped  Polymer  Chain 
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Figure  12.  Isotactic,  Syndiotactic,  and  Tactic  Structures 


35 


ENERGY, 

KCAL/MOiyMONOMER 


Figure  15.  Crystal  Undergoing  Stress 
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.  Force  Constant  Versus  Elongation  of  b-Dimension 


Table  1.  Molecular  Mechanics  Bond-Length  Parameters' 


Bond 

1 

A 

mdyn  /  A 

C-C 

1.523 

4.400 

C-H 

1.113 

4.600 

C-F 

1.392 

6.100 

“Reference  27. 
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Table  2.  Molecular  Mechanics  Bond- Angle  Parameters® 


Bond  Angle 

eo 

\ 

A 

mdyn  A/rad 

C-C-C 

109.5 

0.450 

C-C-H 

109.4 

0.360 

H-C-H 

109.5 

0.320 

F-C-F 

107.1 

1.070 

F-C-H 

110.5 

0.490 

C-C-F 

109.0 

0.650 

‘Reference  27. 


Table  3.  Molecular  Mechanics  Dihedral  Parameters' 


Bond  Angle 

V, 

kcal/mol 

V, 

kcal  /  mol 

V3 

kcal/mol 

C-C-C-C 

0.2000 

0.270 

0.093 

C-C-C-H 

0.000 

0.000 

0.267 

C-C-C-F 

0.000 

-0.086 

0.930 

H-C-C-H 

0.000 

0.000 

0.237 

H-C-C-F 

0.000 

0.000 

0.351 

F-C-C-F 

-0.1000 

-2.000 

0.200 

^References  22  and  27. 


42 


Table  4.  Molecular  Mechanics  Nonbonded  Parameters’ 


Atom 

r 

s 

A 

kcal/mol 

C 

1.9000 

0.0440 

H 

1.5000 

0.0470 

F 

1.6500 

0.0780 

^References  22  and  27. 
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Table  5.  Molecular  Parameters  of  Hasegawa  Structures' 


Parameter 

a-Phase 

P-Phase 

1 

c-c 

1.54  A 

1.54  A 

C-F 

1.34  A 

1.34  A 

C-H 

1.09  A 

1.09  A 

e 

C-C-C 

118.5°,  116.5° 

112.5° 

d) 

C-C-C-C 

179°,  45° 

171.6° 

“Reference  35. 


Table  6.  Energies  of  End-Capped  Helical  Chains® 


mol%  VDF 

Polymorph 

Hasegawa 

"Ideal" 

Minimized 

100% 

diW-trans 

5.48 

5.26 

4.86 

tgtg' 

5.14 

5.63 

5.09 

75% 

all-trafis 

7.33 

7.08 

6.58 

tgtg' 

7.26 

7.85 

7.05 

50% 

all-traTVS 

9.19 

8.92 

8.37 

tgtg' 

9.24 

9.92 

8.87 

®Steric  energies  are  in  kcal/ mol/ monomer. 
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Table  7.  Calculated  and  Experimental  Lattice  Dimensions  and  Angles  of 

PVDF  and  P(VDF-TrFE)^ 


mol% 

VDF 

Lattice  dimensions,  A 

Lattice  angles,  deg 

Crystal 

system 

b 

c 

a 

3 

Y 

100 

8.99 

4.88 

2.56 

90 

90 

90 

orthorhombic 

(8.58'’) 

(4.91'’) 

(2.56'’) 

(90'’) 

(90^) 

(90'’) 

75 

9.08 

5.11 

(2.56“) 

90.01 

89.9 

95.2 

monoclinic 

(8.86“) 

(4.62“, 

(90“) 

m 

(90“) 

5.12^) 

50 

9.10 

5.27 

2.56 

86.6 

87.5 

96.3 

monoclinic 

(9.12“) 

(5.25) 

(2.55) 

(93) 

“Experimental  values^  where  available,  are  shown  in  parenthesis. 

‘’Reference  34. 

‘'H.  Ohigashi  and  K.  Koga,  "  Japanese  Journal  of  Applied  Physics,  vol.  21,  1982,  p.  L455. 

^78  mol%  P(VDF-TrFE);  A.  J.  Lovinger,  T.  Furukawa,  G.  T.  Davis,  and  M.  G. 
Broadhurst,  "Crystallographic  Changes  Characterizing  the  Curie  Transition  in  Three 
Ferroelectric  Copolymers  of  Vinylidene  Fluoride  and  Trifluoroethylene:  1.  As- 
crystallized  Samples,"  Polymer  Papers,  vol.  24,  1983,  p.  1225. 

“55  mol%  P(VDF-TrFE),  low-temperature  phase;  Tashiro  K.,  Takano,  K.,  Kobayashi, 
M.,  Chatani,  Y.,  and  Tadokoro,  H.  Ferroelectrics  1984,  57,  297. 
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Table  8.  Calculated  Crystal  Packing  Energies® 


Hasegawa 

Ideal 

mol% 

VDF 

Polymorph 

vdW 

Coul.** 

Total 

Energy 

vdW 

Coul. 

Total 

Energy 

100 

all-trans^ 

-4.77 

-2.41 

-7.19 

-5.04 

-2.47 

-7.51 

tgtg'  " 

-4.21 

-0.11 

-4.32 

-3.82 

0.36 

-3.47 

75 

all-fransf 

-4.73 

-1.78 

-6.44 

-4.85 

-1.75 

-6.60 

all-frflns“ 

-4.62 

-1.72 

-6.34 

-4.83 

-1.75 

-6.58 

w  " 

-4.33 

-0.47 

-4.80 

-3.92 

-0.12 

-4.04 

50 

all-trans^ 

-4.63 

-1.07 

-5.70 

-4.85 

-1.10 

-5.96 

all-trans’^ 

-4.59 

-1.04 

-5.64 

-4.82 

-1.14 

-5.94 

tgtg'  " 

-4.28 

0.88 

-3.40 

-3.67 

0.89 

-2.78 

“Energies  are  in  kcal/ mol/ monomer.  Crystal  packing  simulations  on  a-  and  3-phase 
super  lattice  cells  were  performed  with  CERIUS  Version  3.1  Software  (Molecular 
Simulations)  for  the  Silicon  Graphics  4D/35  Personal  Iris  workstation. 

'’Coulombic  interactions  were  calculated  using  the  Ewald  summation  technique. 

The  symmetry  was  assumed  to  be  orthorhombic  for  planar  zigzag  ("ideal")  and 
alternatively  deflected  (Hasegawa)  structures.  During  crystal  packing,  the  a  and  b 
unit  cell  dimensions  and  the  setting  angles  were  allowed  to  relax.  Symmetry  remains 
orthorhombic  during  crystal  packing  when  a,  P,  and  y  cell  angles  are  allowed  to 
relax. 

Structure  of  tgtg'  unit  cell  was  assumed  to  be  monoclinic.  Crystal  packing  energies 
were  determined  by  relaxing  a  and  b  unit  cell  dimensions  and  the  setting  angles. 

“Monoclinic  symmetry,  a  and  b  dimensions,  a,  p,  and  y  cell  angles,  and  the  setting 
angles  were  allowed  to  relax. 
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Table  9.  Calculated  Elastic  Compliances  and  d  Constants  of 
PVDF  and  P(VDF-TrFE)  (3  Crystals 


Polymer 

mVN 

(xl0-“) 

Sii'/ 

m“/N 

(xlO") 

S31'/ 

mVN 

(xlO") 

T3, 

N/m- 

(xlO*) 

APs 

mC/m^ 

^33 ' 

pC/N 

PVDF 

2.59 

(8.06^ 

9.52'’) 

0.2 

-0.037 

3.95 

-0.00132 

-3.3 

(-18.8^ 

-25'’) 

75  mol% 
P(VDF-TrFE) 

3.97 

0.19 

- 

2.46 

-0.00252 

-4.0 

50  mol% 
P(VDF-TrFE) 

5.38 

0.2 

-0.034 

1.8 

-0.00190 

B 

^Reference  47. 
'’Reference  33. 
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Table  10.  Calculated  and  Experimental  Sjj®'  Coefficients  and  djj  Constants  of 
Semicrystalline  PVDF  and  75  mol%  P(VDF-TrFE) 


Polymer 

Xc 

%“/  mVN 

(xl0-‘°) 

j  da 

dss  / 

pC/N 

*^33 

(exp), 

pC/N 

8.29 

-10.5 

-32” 

6.50 

-24.8 

PVDF 

^^9 

4.72 

-29.9 

0.75 

2.49 

-23.7 

0.9 

1.15 

-13.1 

9.77 

-10.4 

-34.1  ” 

75  mol%  P(VDF- 

^^9 

7.69 

-24.5 

TrFE) 

^^9 

5.60 

-29.7 

(no  local-field 

0.75 

3.0  ^ 

-23.9 

correction) 

0.9 

1.44 

-13.7 

mm 

9.77 

-12.0 

-34.1” 

75  mol%  P(VDF- 

^^9 

7.69 

-28.5 

TrFE) 

^^9 

5.60 

-34.5 

(local-field 

0.75 

3.0 

-27.8 

correction)*^ 

0.9 

1.44 

-16.0 

=  0.95  from  reference  48. 

‘’Reference  49. 

‘'Reference  16. 

^Correction  estimated  from  Figure  3  in  reference  43. 
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